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results in graphical form on the performance of the receivers 
under different assumptions of wave shapes and p.d.f. on the 
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ET INTRODUCTION 


Classical noncoherent signal detection is generally under- 
stood to mean the detection of a sine wave with random phase 
or time delay in additive white Gaussian noise (WGN). However, 
the extension of this problem to a general noncoherent problem 
involving the detection of an arbitrary periodic signal with 
random time delay in additive WGN, SO ved little attention. 
A case of particular practical interest involves the detection 
of a baseband square wave with random time delay in additive WGN. 

In digital communication systems, the message to be trans- 
mitted is encoded into a sequence of binary digits. Typically, 
these digits represented by the logical states 'l' and '0' are 
transmitted by sending a suitably chosen set of pulses which 
are distorted during both transmission and reception. The 
effect of this distortion is that at the receiver, it is no 
longer possible to determine exactly which waveform was actually 
transmitted. We can model the transmission as well as the 
distortion introduced in the receiver itself, as random noise. 
In all the analysis carried out, it will be assumed that the 
random noise can be modeled as additive, white Gaussian. The 
problem then is one of deciding, on the basis of noisy obser- 
vations, whether the transmitted waveform corresponds to a 
gical '1' or a logical '0'. 

fe ths thesis, the analysis and design of noncoherent 


receivers that optimally (and in some cases suboptimally) 


db 


detect an arbitrarily shaped periodic waveform with random 
time delay, is carried out. The random time delay is assumed 
to have some known probability density function (p.d.f.). 

Two types of p.d.f.'s are assumed for the random time delay \. 
Namely, a uniform p.d.f. and a (parameter varying) non-uniform 
p.d.f. have been considered, and their effect on the receiver 
probability of error (Po) has been studied. 

The design of optimum (in minimum Ea sense) receivers is 
based on the likelihood ratio and the assumption of low signal- 
to-noise ratio (SNR) operations. The design of subopt imum 
receivers is based on heuristic approaches that intuitively 
yield reasonably good performance. 

Chapter II of this thesis presents statistical communication 
theoretic principles that have been applied to the analysis, 
design, and performance evaluation of noncoherent receivers. 
Approximations are made that lead to a decision rule or equiva- 
lently a receiver structure. A justification and meaning of 
such approximations is presented in Appendix B. It must be 
pointed out, that at high SNR, most receivers designed strictly 
from heuristic considerations perform adequately. At low SNR 
however, the receiver design must be optimum so as to not 
further degrade marginal operating conditions. Chapter II 
addresses this issue. 

Chapter III is devoted to the analysis and design of sub- 
Optimum receivers where simple (heuristic) techniques have 


been applied to the design process. Performance evaluations 


EZ 


have been carried out, and wherever possible, comparisons are 
made with optimum systems operating in similar environments. 

The results of Chapters II and III involving receiver Po 
are analyzed and interpreted via the use of tables and graphs 
in Chapter IV. 

The fifth and final chapter presents some general conclu- 
sions to be derived from the work carried out in preparation 
of this thesis, and some suggestions PO reer e analysis in 


this general topic area are given. 


ES 


II.  NONCOHERENT RECEIVER ANACYSIS 


A, < BASIC CONCEET 

Detection of arbitrary periodic signals with random time 
delay is a problem of significant importance in radar and 
communications, which has not received a great deal of attention 
in the past. The most relevant documentation related to this 
problem appears in the radar literature, where the problem is 
formulated as noncoherent detection of a radio frequency sine 
wave burst in the presence of noise [Ref. 1]. 

In communication systems, signals are transmitted to a 
receiver via some medium connecting the transmitter to the 
receiver called the channel. Transmitted signals undergo 
distortion in the channel as well as inthe receiver itself. In 
many cases this distortion is caused by physical processes 
which, because of their complexity must be modeled via statis- 
tical means, i.e., random variables and/or processes. As 
previously pointed out in the receiver itself, noise is 
unavoidably added to the signals to cause further distortion 
and uncertainty. It is often found in practice, however, 
that the uncertainty created by the noise in a receiver results 
in the introduction of signal uncertainties, modeled as random 
processes having known waveshapes but random parameters such 
as random amplitudes, frequencies and/or phases. The problem 
of detecting signals of random time delay has been partially 


considered in Ref. 2. 
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B. PROBABILITY DENSITY FUNCTION OF A 

The type of random signal parameter dealt with in this 
thesis is the phase or time delay of an arbitrary periodic 
waveform of period T. Generally, the period T of the signal 
is very much shorter than the time duration To of the signal 
and therefore it is extremely difficult to predict the time 
delay at the receiver of such a signal. In such cases it is 
reasonable to model the time delay "A, as random variable 
having some density function that reflects the degree of 


uncertainty that exists about \. 


A useful p.d.f. for the time delay [Ref. 3] is given by 
2TA 
mcos- 
_ e 
f£,(A;m) = Im AS [7/2] 2 


O < m<o 


where the independent parameter 'm' determines the degree of 
uncertainty about the time delay A, and I ¿(m) is a modified 
Bessel function of the first kind of order zero. 

For m = 0, we see that Fa (A) - u Foro |IT/2|, so that 
A is uniformly distributed over the period T. As m increases, 
the uncertainty of the random variable à decreases and the 
p.d.f. approaches the shape of a normal density. As m > q, 
£ (A) tends to an impulse so that the time delay uncertainty 
tends to zero and the modeled random waveform approaches that 
of a deterministic or completely known signal. The shape of 
mee p.d.f. as a function of m is shown in Fig. 2.1. 


The derivation of the mean and variance of this density func- 


tion is presented in Appendix A. 


T5 


PDF OF TIME DELAY 


A 





Figure 2.1 PDF of Time Delay 
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ee COMPUTATION OF THE LIKELIHOOD RATIO 

In this section, principles of statistical communication 
theory are applied to the derivation of a decision rule which 
will lead to the design of a receiver that optimally detects (with 
minimum probability of error, Po? a periodic signal with random 
time delay, in the presence of additive white Gaussian noise (WGN) 
under the assumption of low signal to noise ratio conditions. 

We start the problem by considering two hypotheses Hy and 
H, such that during the observation interval (0,T.), under 


0 


hypothesis H we assume that the periodic signal with unknown 


dt 
time delay is present and under the second hypothesis, i.e., 
Hg” there is no signal present. In both cases, the effect of 


the WGN must be considered. Thus symbolicallv, 


Hy: E A yy) 
OE tee le (22) 
— — “o 
Ho: RE) Ate) 
where 
r(t) = the signal at the front end of the receiver; 
v(t) = the T-periodic deterministic signal with 


À a random variable modeling the unknown 
time delay of the waveform; and 


n(t) = a sample function of a white Gaussian process 
with zero mean and two-sided power spectral 
density level No/2 watts/Hz. 


In order to satisfy the minimum probability of error cri- 


terion, the optimum decision rule is obtained by comparing the 


EJ 


likelihood ratio (LR), with a threshold n and decide in favor of 
hypothesis Hj if LR is greater than n, or in favor of Ho 8 
LR is less than n. 

Due to the random variable A, v(t-A) itself is a random 
process. If we condition on A, i.e€., we assume some fixed 
value of A,then v(t-A’A) is completely known. However, the 
observedssıon sr er conditioned on A, because of the noise 
n(t), is a Gaussian random process. In other words condition- 
ing on A, under either hypothesis, r(t) is Gaussian so that 
the likelihood functions can be expressed in closed form and 


are given by [Ref. 4] 


L 


d a2 
WS) Set a (2 
0 o 0 


and 
T 


O 
L, (EA) = F exp {- + f [x (t)-v(t-2’) ] “dt} (2249) 
o 0 


The constant F is of no consequence here as it will 
cancel once the LR, namely L, (x) /L) (xr) is evaluated. Observe 


that 


i) = _J Lj (x |A)£, (A) dA (2.9 
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The likelihood ratio test is therefore 





A 
LAA) 1 
l = > 2.6) 
L. (r) < 9 Ge 
Eege 
0 
where 
n = o la (27) 
PtH} Coy T Ciu 
and 
PtH, } i =. NO (273) 
are the prior probabilities of occurrence of hypotheses H., 
and 
SE HE 1 (2.9) 


are the costs associated with making decisions about hypotheses 


H.. Fora minimum probability of error receiver 
Je EEN 
= 21:0) 
ij | 
0 1 = j 


Taking advantage of some simplifications in writing the LR, it 


is simple to show that the LR becomes 
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I 





L, (x) 2 > O 
nur ” „| RA E (E) (E) dt) 
T Hl 
1 San a 
expi- — f  vó(t-AJatdE, (AJdA Zn (2.11) 
o 0 
a 


Since v(t) is T-periodic, it can be expressed in terms of an 


exponential Fourier series, namely 


2T 


= Jr 
vit) = E E (2210) 
k=-o 
where 
j T/2 -jkt 
SEKR, Í v(t) e def for IS (2) 


Using the Fourier series expansion, we have 


A 5» «@ -j 2 (k+2)) 26 +2) t 
| vitt)d = ) SE e dt 
0 k=-=w L=-0 0 
(2,129 
However note that 
T JÉT (k+2)t | To if £ = -k 
f e se = (22058 
0 | 


0 otherwise 
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provided that TE = nT, where n is an integer. This result is 
also valid even if n is not an integer, provided that n >> 1. 


Thus applying this result to Eq. (2.14) results in 








T 
f ` vÍ(t-)dt = T 5 Ver = y |V E E (2.16) 
one k k SSES k ` 
0 k=-% k=-» 
Observe that E is the energy of thessignal y(t) Tor 0T rn To: 
Furthermore, 
To co são To jkt 
Wee (t)v(t-\)dt = ne r(t)e dt 
0 k=-=00 0 
= VR e (Zon) 
k k 
k=-» 
where 
I: i 
RM | rít)e dt (2218) 
= 0 
So the test of Equation (2.6) becomes 
.277K ei 
oo 2 00 SE 5 À S 
f El ) V,R,e JE, QVGA Con exptE/N_) (2.19) 
00 O k=-» Ho 


SE 


D. TWO TERM APPROXIMATION 
In order to proceed with the analysis so as to derive a 
Signal processing algorithm from Equation (27109) melada 


assume that 


a Ze 
— 1 V,R, e E (2.0200 


In Appendix B, we will show under what conditions this 
assumption is valid. It turns out that the assumption of 
Equation (2.20) is essentially equivalent to an assumption of 
low SNR, i.e., EN, ee 


Since 
es = I +x if nee i 222205) 
provided Equation (2.20) is valid, it às possible to approx 


mate the exponential term of Equation (2.19) by two terms only 


so that the test can be written as 


.2TK EN 
(a RR F V, R u e (A) dd ; TE AN: (2/5208 
EM US A ee o 
=O O k=-» H 
0. 
or equivalenti, 
27k gei 
> V, R r aa E E (A) dr l | (E/N 4-1] E (2.23 
k=-» KK we N < EN e O = Dx : 
Ho 


The test of Equation (2.23) is valid provided Equation 
2720) holds. 

Observe that in order to limit the approximation error to 
ein Equation (2.21), |x| < 0.39. 


In the remaining analysis we assume equal prior probabili- 


À BE E . = e 
ties, namely DATE = — P{H,} and cost assignment or = I Ee 
so the threshold of Equation (2.23) becomes 

N. E/N, 
De le > di (222) 


In order to further analyze the test of Equation (2.23), we 


consider the following two cases: 


12, Uniform p.d.f.- On À 


SE Nenuas Gorm pedi to on ih 


mee Uniform p.d.f. 


For this case, the p.d.f. of A is 


Ea 
£,0) = F E TZ (2.25) 
and the test becomes 
H 
1 1/2 ELA yt 
E VR, f e CO A (2.26) 
k==00 -T/2 
"o 
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However, since 


1/2 27k, pa O 
e di = | (24255 
0 otherwise 


A 
fp | 2.28 
6-0 < Nx (¿alan 
e 
where Vo is just the d.c. component of v(t). Note thar far 


k = 0, Equation (2.18) becomes 


DI 
| 


o P E le (2.29) 
so that the test can be expressed as 
Dad Nu (2.30) 


and the receiver structure takes on the simple form shown in 


pre eae 


Decision 





Figure 2.2 Nonchoherent Receiver (Uniform p.d.f., 
2 Term Approximation) 
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It is clear from Equation (2.30) that this simple receiver 
can only operate provided Vo * O. If Vo is indeed zero, some 
other approach must be used for the purpose of deciding on 
which hypothesis is true upon reception of r(t). In the sequel 
we discuss other detection approaches which will work even if 
Vo = 0. 


The performance of this receiver is evaluated as follows. 


Let 


Zn o (ie) ke Oel) 


Conditioning on A, we see that £ is a conditional Gaussian 


random variable, where 


A 

E{2|H)} O n, 252) 

D 

O 

AR A) = Vo o! a) ete 
a za To e 
= Vo -l Y° = dt 
k=-© 0 
A 
— = 283 
Vor my ( ) 
Furthermore 
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2 
Var(2£|H, ,A) Ef [(2]H, +4) - El2/H,,1)1%) 


T T 


O O 
2 
=. ECV > [v(t-A) +n(t)]dt - Vo o! v(t-A)dt]O) 
T 
2 E GE 
= vá j J E{n(t)n(t)}dt dt = ur ma 
= varítlH,) Ê o? (2.34) 


Equation (2.33) shows that the conditional mean my is indepen- 
dent of A and Equation (2.34) also shows that the conditional 


variance of 2 is independent of A. 








Hence 
2 
FOG) Hee i eee expi —— -~ (2:35) 
d l L l v2TO = 
Q L 
Also, from Equations (2.52) mance.) 
g2 
EE CES = exp{- —s5) (2.36% 
L 0 5 2 
TO o 20, 


The p.d.f.'s of Equations (2.35) and) (2756)) are plotted ni 


Ragen 22.5% 
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T POLJE) 


P(2]H,) 





ES. 2.3 Probability Density Function of 2 


ASE cquations (2.35) and (2.36) and Fig. 2.3, it is clear 
that the two p.d.f.'s are completely identical with different 
means. In order to have minimum probability of error, the 
threshold n should be the point where the two p.d.f.'s inter- 


sect. That is, 





A = 5 (237) 
This threshold differs from that given by Equation (2.24). 
We must however keep in mind that the results were obtained by 


approximations on the LR test of Equation (2.19). 


Receiver probability of error Poe is thus given by 


e EE + E (2.38) 
where 


27 


Pr = Probability of false alarm; 


E 


M Probability of miss; and 


P{H,}, ji =2071. are the prios probabi litiesa 


Since we have assumed equal prior probabilities, Po becomes 




















VOS 1 8/29 1 * vu (mi äi 
ER = 5 e Elk > f E al 
N + V2T0, == 00 (210, 
en do 
co 2 al Q 2 
Br! f A /2ax +4 E f u SEN /2 
N/T» 427 =% /27 
1 Nx n Im ; 
= 5 lerfc+ (—) + erf, ( ) ] (LAN) 
O O 
Q Q 
From Equations (2:24), (22533) ano Ms velas: 
E/N E 2 
já e E = L = Mo - 1 me 
E = CI ] + erf, I - N 1) (2.40) 
2 2 O 
ES N 
N N 
O O 


Observe that if v(t) has a very strong d.c. component in 
comparison to its harmonics, then VAT /N, 1S approximately 
equal to the SNR defined as E/N. - Then Equation (2.40) becomes 
a function of SNR only. Otherwise, 1t Ds a function of bowed 


the SNR and the ver Freezer er 
0 OTO 
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PrZZNonunIrorm pr.d.r. 


For this case, the p.d.f. of A as given by Equation 


(2.1) is 


mcos 
e 


ren) 


Ie is known that [Ref. 5] 


oo 


BROS - Ig(m) +2 J I (m) cospo 
= D 
p=1 
= 1 I_(m)cos pé (2.41) 
See 
p 
since 
Te =e S) 
p ) -p` 


Therefore, Equation (2.1) becomes 


S 21) 
) 1 (m) cos p= 


o 
ELM) = FE, Oma) 





From Equation (2.23) we have the test 


.2TK 
oo — A 


E E UE dA Nt 


aA Vv ER 
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or equivalently 


mE > 
en ee (2243) 
SR 
where 
E -jk 
2TK = Du 
EN ENE dA (2.44) 


Substituting Equation (2.42) into (2.44), we get 





.2Tk 
oi T/2 - J——— 
2Tk 1 2TA T 
EL HA) = == — ) I (m) | COS p = dA 
Ne De TI, (m) SE es p 1/2 T 


7 E T/2 Tp) 
2TI_(m) 1 Tee f e 
0 p=-o P -T/2 


dA 


T/2 57 (p+k) A 
8 f e dA] (2046) 
-T/2 


Note that both integrals in Equation (2.45) are zero unless 
p = tk. For p= k, first integral in brackets is T, and the 
second integral becomes zero. For p = -k, the first integral 


is zero but the second integral is T. 
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Therefore 








I m) 
Zaks = k 
En (A) - I, m) e (2.46) 
since 
I y. ) = Ty e 
Hence the test of Equation (2.43) becomes 
T > 
I pla WKD a SC? 
0 Kam H 
0 
The test can be simplified further by letting 
“= e 7 V, Ry, I, (m) 
0 k=-» 
T T ZN 
- d Pt Am E 
= ——-{V,I, (m) ee m r e 
E SO y kl "TP o 
E TA jet 
+ ) VI (m f (ae dt) 
al 
-k=1 
T 
O 
= mW I, (m) u ae 
Al 
> 2 27k 
+ 2 A |v_ II, (m) o! cos o dtl (2.48) 
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where 


tja, 
Vu AAN (2.49) 


Finally the test becomes 


T . T an 
Í f ee (EE +. 1 
vo a Se l e (m) r(t)cos(==t +a, > ie 
0 Ore k=. © 0 
0 
(2.50) 


and the receiver may be implemented as shown in Fig. 2.4. 


2|v, Ir, (m 


In (m) 


2|v,|r,(m) 


em — 
I) m 





Figure 2.4 Noncoherent Receiver (Nonuniform p.d.f., 
2 Term Approximation) 
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Notice that this receiver can operate on signals both 
with and without d.c. component. When a signal has zero d.c. 
term, the receiver of Fig. 2.4 can be simplified by eliminating 
the uppermost branch. 

The performance of this receiver is obtained by using 
the fact that conditioned on A, under either hypothesis, g 


mona Gaussian random variable with 


Elg|H,) = 0 | (rat) 
and 
T Rene 
o é S e 
Elg|H à} = a A. V, 2), (m) o v(t-A)e dt} 


d 00 oo 
T tay KE BN lk f e dÈ 


(252) 


The integral of Equation (2.52) has been analyzed in Equation 


(2.15). Therefore 


h (A) (2253) 





T 
Elg|H,,A) = ) |V,|"I, (mje 


Furthermore, 
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T SE 


© O 
dk 
Var(g|H,,A) = egen Wd ™ 2 IV (ET e dt 


xo i SE 
1 TE > 
E Tm) o! vít-A)e dt) } 


= To j@uk, 
E .W,1,(m) TI n(te 7 at)*} 
= —00 O 


E, 


T 
CO co 


1 O 
= y l V V,I (m)I, (m) ff E{n(t) niga 
Fal) Jee 0 








No : : d Stir 
> ae o Lo mom) e dt 
I, (m) 2 Kal gato E LK S 0 





EES 
= oe Vee 
213 (m) k==w K 8 


A 2 
NAO O o AS Io (20505 


With p = P{H,} and 1-p = P(H,), we have an expression for the 


conditional Po given by 
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Nx 








P (A) = p E a (I-p) J E(g|H À) dg 
* 
oy =g*/20% Me =(2-h(0)) 2/207 
= p | —e “dg + (1-p) f e “ag 
EE vio =% /270 
g g 
Nx NESTA (A) 
= perfect (petit  ) Eos) 
g | g 


where n,, h(A) and SÉ are given by Equations (2.24), (2.53) 
and (2.54), respectively. 
The actual receiver Po is obtained by integrating 


Pa (A) ever the p.diFf. of A, resulting in 


a o AO 
P, = p Pen + (1-p) J og S (2.56) 
Observe that 
SE E/N 1922 
O O 2 
Nx q (e 2l) 
a NT co 
g 22 ) Ev, 1 1 (an) 
217 (m) k=-0 
0 
E/N, e 
E = E — r E 
2E 22 2 2 
EL t (m) To (m Bak | 
Also 
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2 SE 








T © J—— 2 
52 vêr (me >, 1/2 
Ii,.(m) ko 
EONA 0 
g NT co 
g O QO y |V [212 (m) 
Oo a e 
0 
= -2 NK z 
= rina me a A 
=- ORO TES A A X (2.58) 
E Dim Tv? = | 
NE a a: 0 E 


E. THREE TERM APPROXIMATION 

whenever Equation (2.20) is not satisfied to the degree 
that Equation (2.21) is valid, it is possible to proceed with 
the analysis beyond Equation (2.19), if in Equation (2.19) we 
expand the exponential term and approximate it with three 


terms. Since 


exp{x} 2 1+x+--, if x <<] (20580 


in order to have the approximation error not exceed 10%, we 
must have x < 0.79. Using the chree cermappros metendo 
Equation (2.59) and applying it to Equation (2.19), the test 


then becomes 








A a _.2rk S _42mk) = 
2 T" 1,2 EE > E/N 
J [1 KS EA RS ty AR “je, (ar me 
o E H 
0 
(2.60) 


or equivalently 


‚2Tk ah 


oo co o NM co co ai 
2 JT il 2 ql» 
= E De f £, Je E f nz RH V,R,€ ) 
O k=-0w =% -%0 O k=-o 
a 
> 
E (AJA É Y (2861) 
Zo 
where 
E/N 
y = (ne erh, 262) 


Here again, two cases must be considered: 
ime Uniform p.d.f. on A 
Em Nonuniform p.d.f. on À 
Uniform p.d.f. 
For this case the first term on the left hand side (L.H.S.) 


of the inequality in Equation (2.61) has already been analyzed 


Mec uations (2.26) through (2.28). This first term is just 
VoRo: The second term of Equation (2.61), i.e., 
E S 421K, 
1 2 
> SE I vre * )% 009 


becomes (with the aid of Equation (2.15)) 


Be 27 
œ 00 -j—(K+2) A = 
l 4 T 2 
35) Y V,R,V,R, f e da == > |V 
2T y? ee to Es LA N2 k 


Ee 


The test therefore becomes 


R CO 
1 2 
N a pv IR 
O 





yA (2.604) 


k=-» 


or equivalently (using Equations (12.15 ane eco 


T T 





Vo O Vo O 2 
A oe 
o 0 o 0 
H 
1 
o IM. I 2 
ke E É ° r(t) cos m ae j yF (2005) 
LON N T < 
k=1 D o 
Ho 


This test leads to the receiver structure shown in Fig. 2.5. 
Note that this receiver utilizes both the d.c. and non- 

d.c. components of the periodic signal. However, in practice 

there are many cases where the d.c. component of the signal is zero. 


For those cases, the test reduces to 


e (2.66) 


and the corresponding receiver can be implemented by simply 
eliminating the upper most branch of the receiver in Fig. 2.5. 
2. | Nonuniform pica 
Our starting point here is Equation (2.11), repeated 


here for convenience, 
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ect) 
cos Ze 
e 
E hold for 
o comparison 


e e 
e o 
e © 


Fig. 2.5 Noncoherent Receiver (Uniform p.d.f., 
3 Term Approximation) 
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T T pal 
j 2 2 l = > 
J expl J ritjvitojdtjexpi- — ) vw tt-A)jdere Vor 
=% o 0 OTO y 
0 
where n is given by Equation (2.7). Note that the second 
exponential of the inequality given by Equation (2.11) has 
already been analyzed in Equations (2.14) through (2.16). 
Furthermore, if we expand the first exponential in the above 
inequality and approximate it with three terms as given by 
Equation (2.59), then the test of Equation (2.11) becomes 
T 
| Y fo rtt)vit-adt £, (AJdA 
N A 
. De el 
+> f tA fo ríe)vte-ajatlóE, (A)dA 2. y (2.67) 
A < 
—00 o 0 | H 
0 
where y is given by Equation (262, 
Interchanging the order of integration in Equation 
(2.67), we obtain 
T 
2 S ` 
wf een f/f (ta) £, A) dadt 
o 0 =% 
T E St 
2 3 > 
+ > ite) Í v(t-A)vitT-A)E (A)dA dtdt Cy (2.68) 
No 0 =00 Ho 
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Note that v(t) is a deterministic signal but as mentioned 


earlier, due to the r.v. A, 
A 
Y (E) MASA E) (2200) 


itself is a random process with expected value 


oo 


Ely(t)) = ff vitor (da É mi (t) (2.70) 


= 00 


and autocorrelation function 


oo 


E{y(t)y(t)} = f V(t-A) v(T-A) Ey (A) AD 


= CO 


A Belle a) 2.71) 


Substitütion of Equations (2.70) and (2.71) 
in Equation (2.68) leads to the following test that a receiver 


has to perform 


E T To al 
= o! r (t)m, (t)dt + + f le r(t)R, (t,t) Jdt i LO) 
O o 0 Ho 


and the implementation of this receiver is shown in Fig. 2.0. 
Here the autocorrelation function Ry (t,T) can be considered 
as the impulse response of a filter that has to be designed 
for detection of the signal. 


Purcnermore, for the unirorm p.d.f. On A, the expected 


value of y(t) given by Equation (2.70) becomes 
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Decision 





Fig. 2.6 Noncoherent Receiver (Nonuniform p.d.f., 
3 Term Approximation) 


AS = ae 
-T/2 
, Sg 7/2 Ja 
== Ll ve foe dA 
k=-» -T/2 
z Vo (2:18) 


which is a constant, and the autocorrelation function ota 


given by Equation (2. curas GUL to se 
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> E: ¡e ¿e T/2 ÍSE(K+L)A 
ert) = = } À V,V,e e f e dA 
km 2=-0 -T/2 
E y IRE 
SEN E (2.74) 


emiten is a function of (t-t) only. This implies that for the 
fetter p.d.f. on A, the impulse reSponsenof eheytilter turns 
out to be time invariant. That is, it depends only on the 
time difference (t-t). Furthermore, my (t) being a constant 
can be eliminated from the upper branch of the receiver shown 


nro. 2.6. 


F. ANALYSIS OF TWO SPECIFIC SIGNAL WAVESHAPES 

as ipouncediout earlier that the receiver derived using 
a two term approximation on the exponential appearing in Equa- 
tion (2.19), is unable to discriminate between the signal 
present versus the signal absent case for signals having zero 
d.c. component when the random time delay obeys a uniform 
distribution. When the random time delay obeys a nonuniform 
p.d.f., a receiver was derived (see Fig. 2.4) that could 
detect the presence of signals having zero d.c. component. 

Machis section, assuming a nonuniform p.d.f. on A, the 
receiver performance will be analyzed for two different signals 
which are found quite often in practice and that have a zero 


d.c. component. The two signals are 
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1. Sine Wave 
2. Square Wave 
1. Sine Wave 


For the first case let 


were) 


Il 
> 
Q 
O 
LD 

| 
ct 


(2273) 


by 


This periodic signal has only two discrete components. 
In other words, it can be represented in terms of its two 


Fourier coefficients, i.e., 


NJ > 


V = | (2:70) 


0 otherwise 


Thus, 


in) = (22573 


and from Equations (2.57) “aadr (2.5 Meño Eat 


E/N, 
J ae O (2.78) 
Cg i E, 1072 l 
T 
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E/N T m) 
BG L) 
Igo ™) E 
ae (A) O (2.79) 
o 7 I, (m) 1/2 ; 
3 E 
Thus the expression for Po becomes 
E/N, 
Ee ee Ee (ne -1) 
S a 12 
a 
0 
E/N ES (im) 
O 265 l 2TA 
(1-p) 772 nen = nee 
ae o 0 T 
eh mn / EU e EE AA 
0 E 
“Te A NO A 

I ¿(m) O 022.80) 


Computation of Equation (2.80) can be ee, out as 
a function of the SNR E/N.) p and m. This has been carried out, 
and the results are presented graphically in Chapter IV. For 
the sake of computational simplicity,equal prior probabilities 


were assumed, i.e., 


Besultıng in n= l. 


Dr, Square Wave 


For the second case, let v(t) be a periodic square 


wave as shown in Fig. 2.7. 
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Observe 


and 


Eg. 


2.7 Periodic Square Wave 








The Fourier coefficients of v(t) are 

ee 327%: E ¡E 

mi f Ae da Ae dt] 

0 722 
Sin TK 2 
jak =3Tk 2 
>= Ae [Tk (2258718) 
2 


that ve and all the even coefficients are zero, i.e., 


Vox = 0 for all k (2502) 
ee wis 2 2A 2 
= 5 prio = oe MA (2088) 
k odd 
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Note that 


T/2 > 
f ve(t)dt = AT = T ? A 
-T/2 k= =00 
which implies that 
a 2 2 
a zen 
Kan 
Furthermore 
E 2 2 S 2 
IL (mm) = |v_ (1. (m +2 v Im) 
so that 
ee 
7 |w, 121, (me E ze IV | I (m) 
foe x 0 | 
oo 2rk 
+ 2 Ge E 21 „ (m) cos A 
k= 
since 
O E I, (m) =  I_,(m) 
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(2:84) 


(20853) 


(2 86) 


VANA 


Hence Equation (2.57) becomes 


E/N 
(ne 1) 





Nx 
— (2783) 
k=l 


A“1ó (m) 


3 4E/N |) Lv, | “x2 (mn) 


and Equation (2.58) becomes 


E 
N_ 4E/N 
(ne °-1)-(-— EI oT) 
= e E EE (2.89) 
g CO 
4E EE 
— Pv CT (m) 
No VEH k k 


15 (m) AS 


2 
| I, (m) cos 





Thus Equation (2.56) for the receiver probability of erros 


becomes 


E/N, 
P = p erfcx une =) 
2 
4E/N, Lv 1, (m) | 1/2 


2 2 
E A 
q 


E 
N 4E/N © 
2 27k 
(ne °-1)-(——>-5) ( J _ |v, | “1, (m) cosa) 
T I,(mA k=1 
+(1-p) | erf, 
0 4E 
™ 2 ie 
oO k=l 


2 2 
I A 
gm) 


a SE 
| I, (m) | 1/2 


meosZ TA 
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Observe that Equation (2.90) (as well as Equations 


(2.57) and (2.58)) involve summations whose index runs from 
l to infinity. In practice, it is not possible to compute 
infinite sums. However, due to the fact that for all reason- 


able signals v(t), the magnitude of the V, coefficients gets 


k 
smaller as k gets bigger, the infinite sum can be truncated 
without introducing significant computational error. Further- 
eee, terms involving the magnitude squared of these coeffi- 


cients will tend to decay more rapidly allowing early truncation 


of the infinite sums. 


6% DISCUSSION 

In order to carry out the evaluation of Equation (2.90) 
to a reasonable degree of accuracy, the mathematical expression 
given by Equation (2.90) was evaluated by computer for fixed 
SNR values. Computation of Po was carried out for a fixed 
value on the upper limit on the index k. Recomputation of 
P_ was carried out everytime this upper summation limit was 
incremented by 1. Recomputation of Po was stopped after an 
incrementation of the upper limit in the sum terms did not 
yield an appreciably different value for Es in comparison to 
the P value before incrementation. 

The result of these evaluations has been presented in 
Chapter IV in the form of tables and graphs while P has been 
plotted as a function of SNR for different values of the 
parameter m. Recall that m specifies our relative prior 


knowledge about the signal's random time delay. 
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III. ALTERNATIVETSUBORTTMUMFREBEHVERS 


In the previous chapter, optimum detectors based on 
likelihood ratio tests were derived under conditions that are 
equivalent to low SNR conditions for signals having random 
time delay in the presence of additive white Gaussian noise. 
The tests derived gave rise to different receiver structures 
which, whenever possible, were analyzed in terms of their 
error probability performance. Both the receiver structures 
and their performance were seen to be a function of the uncer- 
tainty in the signal time delay or more precisely, a function 
of the p.d.f. of the signal time delay A. 

In this chapter, two suboptimum receivers are analyzr ce 
The approach used here is different in the sense that re- 
ceivers are derived via strictly heuristic means and their 
performances are then analyzed. These receivers are: 

A. N-correlator Receiver 


B. Estimator-correlator Receiver 


A. N-CORRELATOR RECEIVER 

This receiver is shown in Fig. 3.1. The idea is to corre- 
late the received signal r(t) with N delayed reference signals 
generated at the receiver. If a signal is present, presumably 
at least one correlator branch would produce a large output so 


that the sum 2 could be guaranteed to be much larger under 
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Fig. 3.1. N-Correlator Receiver 


these conditions than when the correlators are fed by low 
level noise only. Thus, if 2 is greater than a prescribed 
threshold, then the presence of a signal with a random time 
delay is declared; otherwise only the presence of noise is 
declared. Here we assume for convenience that the observation 
time To is an integral multiple of T, the period of the 
signal v(t). Observe furthermore that we are not required to 
assume low SNR conditions as done in Chapter II. 

In order to analyze this receiver in more detail, observe 


that 


le 


N O aa N 
i PS II E 
n=l 0 nel 2 
where 
T 
A z AT 
eo = E EE (3.2) 


If we consider the two hypotheses, as done in Chapter II, 


namely 
Hy: r(t) = v(t-A) + n(t) 
US “EMS To (3D) 
Ho: r(t). = nE) 


then conditioned on hypothesis Hu, £ is a Gaussian random 


variable (r.v.) with 


Elr Ho) = 0 (3.4) 


since the noise is assumed to be zero mean, and 


T T 
E = ad i P f TESTER USERS 
nl ee A 2 u É N 
T 
N O N 
=. ae ot D = 20) 
SC R| Re od VAR ehe as (3455 
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where 


T 
A S nT mT 
en e Ee (oo) 
are the signal discrete cross correlation terms. Thus, 
N 
ee O (SU) 
n=1 
and 
í ~ E 
Eto |H.} = Ei E GE 
q n=l m=1 Nm 0 
N N N 
8 O A 2 
Ss en nn 
n=1 m=1 


Now conditioned on both H, and à, £ is also a Gaussian 


il 


Eevee with 
o (A) (3.9) 


Etr HAD = f  v(t-A)v(t - )at 


and 
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N 
Be) 
and (lie Ele, [4 AAA | o (3110) 
so that 
N 
E(2/H,,A) Si W a, (A) (32135 
n=1 
and 
El [2-—E(2 |H Ma A} 
De 12 
N "o ls "o ES D 2 
= Ell  (f vít-A)Jvit-Jdt + { n(t)vit-— dt) E 
n=1 0 0 n=1 
T 
N O 
EE ert? 
n=1 0 
zi? 
Ee A = 0 (32128 
2 Bene nm £ 
Thus the conditional p.d. S Or Rare 
1 -24/207 
Ep (LH) = —o (3.13) 
Jonoa 


2 


and 
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a E 
ES LO) GER 
Ii n=l 


£ (2H, 1A) = ———e (3.14) 
Varo; 
with 
E (IH) = E E (HA) E, (AA CORES) 


me detector of Fig. 3.1 performs the test 


A IS 


so that its error probability, Po is given by 


CO 


e 
m A p, (2]H,142 + P(H,)  f p,(2/H,)a2 (3547) 


= OO 


Let p = Pina and 1-p = P(H, ) so that using equations 


(ees) and (3.15), equation (3.17) becomes 








E 2 
-(2- Y a (A)) 

n=l 

2 2 2 

co -9/20 Y 00 20 

PL = p f—e ‘anv(l-p) ft f ——e E 
y v270 -> -=> 270 
A x 

£¡(M)dAJ al RS ta) 
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Interchanging the order of integration in the second 


(double) integral of equation (3.18), we obtain 





E 2 
O 
n=1 
2 
© Y 29 
f tay| f —e É az | aa 
— 00 =% /210 
ie 
N 
EE 
n=l 
O 
co 2 2 
= Dee í E TE a 
= OO — OO vd Am 
N 
a ML Ay 
n=1 
= e EHNEN (3.19) 
O 
em OO 2 
Thus, the receiver probability of error becomes 
N 
A e a ON 
= a E u N a 
Pa = P A + (1-p) = Ee A dA (gay 
The threshold y can be optimized to minimize Po: Such an 


approach unfortunately leads to a complicated integral equation 
that is not easily solvable. Tomor ek rennen ch) E can be 
calculated on the computer for y = 0, and then incrementing 
|y| until an optimum value of y is found for a given SNR. 

It must be pointed out that the receiver of Fig. 3.1 will 


not always yield desirable results. For example, if v(t) is 
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sinusoidal, in the absence of noise, % = O which means 
that with noise present, the receiver is unable to detect 
He presence of the signal. Fortunately for v(t) sinusoidal, 
a better receiver is available and its performance is well 
known [Ref. 4]. 
l. Special Case of Triangular Wave 

Consider now a periodic waveform of triangular pulses 
and its two arbitrarily eee versions as shown in Fig. 
3.2. If we represent one triangular pulse by p(t), then we 


can express v(t) as 
M 
Dee O (sa 


In order to find the cross correlation terms N 
required for the computation of Po lequation (3.19)], 

instead of taking all the pulses in equation (3.21), we 

simply consider the cross correlation of a single pulse of 
period T with a delayed version of itself and then multiply 
the result by the total number of pulses M within the interval 
To’ since v(t) is assumed to be periodic with MT = To 


Gëfter: Cen make computation simpler, Fig. 3.3 shows p(t), 


p(t-a) and p(t-8), where a < B has been assumed. 


Sr 


v(t) 





MALT MT 

v(t-a) 
@ ® e 
E 

0 a T+a ¿T+a (Mel) Tra MT 

v(t-8) 
$ t Oo Ñ 
0 8 T+8 (M-2) T+8 MT 


Fig. 3.2 Periodic Triangular Wave and Two of Its 
Delayed Versions 
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p(t-a) 










AT 
EE 


At+A(T-a) 


A(t-a) 


4 Y (8-8) 
A At+A(T-8) 


A(t-8) 
t 


Fig. 3.3 One Period Restricted Triangular Wave and 
Two of Its Delayed Versions 
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Observe that for a < 8 


T 


O o 
f v(t-a)v(t-68)dt = M] f p[t+(T-a)]p[t+(T-8) ]dt 
0 0 
B IF 
+ f p(t-a)p[t+(T-8)]dt + f p(t-a)p(t-8)dt 
a B 
a 
= M| f [At+A(T-0)] [At+A(T-8) ]át 
0 
B T 
+ | Alt-a) [At+A(T-B)]dt + f A(t-a)A(t-B)dt 
a B 
2 2 2 2 3 
e? Se En E TASA T 
a E 2 * RE 
(ION 
POr Q3 e and ß = N, equation (3.22) becomes 
2 2 
2 30 n nm m m I 
MAT’ (— + = - — + == +=») = A (323) 
N 2N q on? 2N 3 nm 
provided n < m. For n = m or equivalently a = 8, the quantity 
inside the parenthesis in equation (3.23) becomes 3º 
Similarly for a > B, equation (3.22) becomes 
o 2 2 2 2 3 
J v(t-a)v(t-8)at se WU SS agr +20 48 +40) (3.24) 
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Letting a = Ze and ß = = in equation (3.24), we obtain 


= 
5 
a 
3 
pa 


A = MATA 4 +50 45 (3.25) 


N 
Z 
Z 
N 
Z 
N 
N 
Z 
W 


provided n > m. Here again for n = m, the quantity inside 
the parenthesis in equation (3.24) becomes >. 
To proceed further with the analysis, the double 


summation of equation (3.8) is broken into three summations 


meamaccount for n < m, n=m and n > m. Thus 


Kr M y ES 
A = A + A + A (3.26) 
n=1 m=1 dai m=n+1 n=1 N naa PR n=m+1 m=1 ek 


Observe that the L.H.S. of equation (3.26) represents 
the sum of the elements of an N xN matrix. Whereas the middle 
term on the R.H.S. of this equation denotes the sum of the 
diagonal elements of the N xN matrix, the first term is the 
sum of all the elements in the upper triangular part and the 
third term is the sum of all the elements in the lower triangu- 
apart of the N xN matrix. 

In order to evaluate the first term on the R.H.S. of 


equation (3.26), use of equation (3.23) (n < m case) yields 


6l 


Nal N NSI N 2 2 


223 n n nm m m dl 
A = MATT” } [— base - — + -— +5 
DST nn n=l m=nt¢+l 20° N wi on 2N 3 
N-1 N N-1 N 
- eee po PRM 
2 4 u 2N 4 = 
2N n=l m=ən+1 n= 1 mani 


1 Nt ; 1 Nt l ; > 
= = l nm + = m 
N“ n=l’men+L 2N n=l m=n+l 


Ei 
2N 


2 
E Maq A (3.27) 


The second term on the R.H.5S) of equation LS) 


becomes 
N N 
L Ae = MATO ] q = mul (3.26) 
n=1 n=1 
since A SE for n =m 
mn 3 = 
Finally, substitution of equation (8225) Nin > m cas 


in the third term on the He eS 


N N-1 2 2 
23a n nm m m JL 
) A = MATT" [— - 5 +>5 +50 +5] 
n=m+ I m=1 DIE on? 2N Nº 2Nº 2N 3 
DEE Ne 
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Addıng the resmlts el equations (3.27), (3.28) and 


fees) yields 


N N 2 
p 2m3,3N" +1 
A ae ae! (3.30) 


n=l m=1l 


Thus fromiegquations (3.78) and (3.30), the conditional 


variance of 2 becomes 





No 5 y A eee 
2? oi ni DÉI a GE 
N 2 
E O SIM pd 
e (ars) 
since 
E. >» e 
fee ct )ct = M {vi (t)dt = M f (at) at 
0 0 0 
2,3 
_ MAT Å p (3.32) 
3 
ee energy of v(t) for 0 <t < To: 


In order to compute the probability of error Por We 


N 
need to evaluate ) a, (A) given by equations (3.9) and (3.11) 
n=l 
as 
N N a 
WER 
) a (A) = ?) ft  v(t-A)v(t -=7)at Gase) 
n N 
n=1 n=1 0 
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A Fourier series expansion of v(t-A) and v(t gol) 





N 
yields 
N nn ~ 3 KA “jo "o jr +) 
) EE ) ) ) HK e f e dt 
n=1 n=1 k=-wo L=-0o 0 
2 re 
o0 > Gx N jn( LES 
= T ve ee ) e (eur 
=-o n=l 
Observe that for k = +iN, where i is an integer, the 


second summation involving the index n yields N, otherwise it 
yields zero. 


Thus, equation (3.33) can be written as 


N PR > -jkm 
1 SE TN (Ivo + 2 IS 
n=1 k=-» 
a , Ik 
+) gal € 
GN 
= TNV]? +2 J ul" eps Fam (3.35) 
k=l 
since 
2 2 
IV, | = Na 
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In equation (3.6), if we represent v(t I) and v(t I) 


in terms of Fourier series, then 











Sa E Ea A T Gert 
Me an = OU dO YY € f e dt 
n=1 m=1 n=l m=] k=-~© l=- 0 
2 o 
= > N "nt: ZS N jm( =s) 
= To } IV, | } e e 
k=-» i als 
S 2 |sin(kr) É 
= T J ml S- (3.36) 
k=-o sin(—T) 
since 
N no sin(N9/2) _3(N+1)0/2 
2 E ESA ~“ a 


Observe that the R.H.S. of equation (3.36) is always zero 
except for k = +1N where 1 is an integer. Also observe that 
for k = +iN, the R.H.S. of equation (3.36) becomes an indeter- 
minate form. Due to that reason, differentiating sin(k7) and 
ein? w.r.t. k and then taking the limit k > 1N, equation 


(3.36) becomes 


; ` E i ` 
A = TN IV, | 
nel mel MM SNE 
2 SE E 2 
m „L Yon! ) (esse) 


so that equation (3.8) can be written as 


2 O 


Sa E (369 


Equations (3.35) and (3.39) are of general form and can 


be applied to any signal, not just a triangular waveform. 


The Fourier coefficients Hi of the periodic triangular 


wave are 








AT BR 
2 en 
Ke = (3220) 
2 otherwise 
and 
22 
ae k=0 
vl = (3.41) 
22 
AT . 
775 otherwise 
4n k 
so that the use of equations (3.31) and (3.35) yields 
* 
528 A ; n Y (3% 42) 
O 2 2 
É m ee 
O 8 N 
O 
where y* = Y/N, 
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21 
N 22) air“ cos km N 











o (A) N( +2 J um, 
n=l 3 B E 8 ) 4 k=l Amen? a 
ER No 30° +1 AT“ ne 
o aa aa 
Kel Ay e 
wil +4ji-44 (221) 
E 3 N? 6 T T 
0 3N" +1 1 + 1/3N 
since 
00 2 2 
y *5* - E-lo+£ 0< 0 < 27 (3.44) 
Sl k 
and 
N 
Y- ba (A) 
n=1 " E w E 8 
= - 2 RO 
g, B 3N +1. 1/2 o 3N +1 
[ Aeg h 
O 
SES A eZ 
N{l oh, T t (m) 5 
d —— > —Á (SAS) 
l + 1/3N 


Finally, the expression for Por after substituting equations 


ea, (3.42) and (3.44) in equation (3.20), becomes 
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nee ÓN 
BE T 
P = p erfcx ÓN (1-p) 
= y E GN “+1),1/2 -T/2 
N omen 
O 
er ft =e ee 
¡E LAN +1),1/2 o 3N%+1 
N 8 
3 Te 
3N = a I l. 
dÀ (3.46) 
(3N? +1) 


As mentioned earlier, the receiver of Fig. 3.1 has 
been designed via ad-hoc methods, so that the optimum threshold 
setting is not readily obtainable. The difficulty in applying 


an analytical approach to the problem of finding the optimum 











dP 
threshold arises in the solution of u. = O for y. Since 
N 
(y- Y a, (a))Í 
n=1 
aP. i -22/204 T/2 i 20% 
a PS te) E dÀ 
V2T0, -7/2 d Zo, 
(3.47) 
dP 
then, solutionter a = 0 is equivalent to solving 
I ap 001? 
Y UA) [ GA) | 
del _ nel S 
WZ os 20% 
SIE e dk = Tê) (3.48) 
-T/2 p 
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Clearly, the computer could be used to search for a value of 
y that satisfies the equality prescribed in equation (3.48). 
The approach that has been used in trying to obtain 
the optimum threshold involves using the computer to evaluate 
equation (3.46) as a function of y for fixed SNR. Thus the 
value of y for which Pe is minimum yields the optimum thres- 
hold for that particular value of SNR. This procedure is 
e peated for different values Gt SNR.» the resulting values 
echen used to plot Po VS. SNR. These curves have been 
plotted using optimum threshold settings and the results have 
been summarized and discussed in Chapter IV. (Observe that 
the threshold y is SNR dependent. This situation is typical 


of noncoherent signal detection problems.) 


DEE STIMATOR=CORRELATOR RECEIVER 

The second suboptimum receiver being proposed is the re- 
ceiver structure shown in Fig. 3.4. It is basically an 
estimator-correlator receiver in which first a coarse estimate 
of the time delay in v(t) 1s made and then this information 
is used to process the delayed signal r(t) coherently with 
an estimated reference. 

In this receiver the largest output of the N correlator 
branches is used as an estimate of the signal phase by the 
logic choice. This estimate is used to provide an accurate 
local reference with which to perform the correlation detection 
Operation. The output & of the main correlator is used for 


threshold comparison in order to make the decisions. 
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Thus 


T 
O 


2 = fo r(t)v(t-a)dt (3.49) 


NO 
GA 


I Tar e 
Nº N'N' 


value for 0 SES: d 


where a = +: T and «a may take on only one Such 


= In order to determine receiver perfor- 
mance, we note that 2 is a conditionally Gaussian random 


variable with 


E(2|H,) = 0 (32509 
T 
2 > 
vVarí2]H,) = Ell f  n(t)v(t-ajdt]*) 
0 
N To 2 
= > fouvílt-a)dt 
2 
0 
No oo co ¡+ a = jfi (k+) t 
= om VV oe f e ae 


LG 

E 

8 

Z 
ID 
N 





E 
= > E O El 


where Vi are the exponential Fourier series coefficients of 


v(t), and E is the energy of v(t) for 0 < t < T. 


Dar 


Also 


E{2|H,,A} 


JE 
O 


f Em Oy GE montei 
0 


SE (kara) = Tx+2) e 


) E Ä f e de 


= == 00 


2 = 2 2 A 
Zaun +2 } IV, | cos Tk (A-a) ] = 


It can be shown without a great deal of difficulty that 


2 


Var(2|H,,A) A 


Therefore the conditional p.d.f.'s of 2 are given by 








1 -27/205 
(2|H,) = e 
Q /270 
2 
1 -[2-m (A) 17/204 
(2|H,,) = e 
(210, 


(3082) 


(ISSN 


(SAM 


(ION 


SO that 








co 1 -[2-m(à) 1/20% 
o = f e E) da (ere) 
L 1 = A 
=% Y2TO 
Q 
Thus the expression for Po is given by 
>. Pi E E lash se) 
>,  -22/20? e -(2-m00)/20% 
= p — e dg + (1-9) | [—— 
Y 270, -0 >=% V2TO, 
x £, (A) dade (2257) 
where 
AS and (l-p) = P(H,) 


as before. 


After a change of variables, equation (3.57) can be 


Written as 


l _-x”/2 Be 
P, = P f —e +(1-p) f/f f =e dxf, (A)dA 
E Y/O, /2T -o —o /27 A 
= perfes(L) + (1-p) f 8,0) er aD ay (3.58) 
To Bes A To 
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Since the receiver was not designed from likelihood ratio 
test principles, the threshold y must be set for minimum Po: 


The approach used to find the optimum threshold involves 











dP 
solving aa 0 for y as before. 
Since 
* 
+ = Lo = E ere EE (325 
Q /N_E/2 VE/2N_- = 
O O 
and 
|V pe +2) IV E ce 
m(A) E k=1 É 3 
EAN m e aieeaa ( 3260) 
Si ° 2 S 2 
aal + 2 Do lvl 
we have 
[vo |? +2) Iv, | ºcostTE(a-a) 
= * _ = 
ul a SB V2E/N_ E, AE (3,68) 
2 VE/2N vl? = Iv, 1? 
Thus equation (3.56) becomes 
mcos HZ 
E E T 
P= p erfcal— — ) tr (pi or 
q /E/2N ec 
O -T/2 
|V Ee N cos EME (7-a) 
N o 0 er T 
x erfx —— - V2E/N ~ dr 
E/2N 2 2 
VE/EN, vol? +2 Flv 
L=1 
(3.625 
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As before, determination of the optimum threshold by solution 
AP 
of — = 0 involves solving an equation of a form similar to 


dy 
that of equation (3.47). 

In general, solutions in closed form do not appear tracta- 
ble. A search for an optimum y may be performed or solutions 


for specific cases may be possible. This is illustrated by 


considering two special cases in which v(t) is a 


ay sine wave 
2. square wave. 
l. Sine Wave 


As a first special case, let 


V(t) = Acos en, 
T 
Bon 27 
e En 
so that 
A 
= .64 
Vi | (3264) 
0 otherwise 
and 
© 2 2 
2; vi = 2% = 4 CE 05) 
k=l 


Thus equation (3.60) becomes 


A 27 
E en cos 21 
= vV2E/N — = V2E/N_ cos zn, (3 616) 


Si o Al /2 














and for p = (l-p) = > equation (3.58) becomes 
qu E 
E = 2 SO a ote TI_ (m) f = 
VE /2N, 0 -T/2 
* 
x erf = v2E/N qa (3160) 
== O T 
VE/2N | 
O 
As mentioned earlier in Chapter II, m = 0 corresponds 
to- a- uniform p.d. ona ee e so that the error 


probability expression for m = 0 becomes 


| x T 2 + 
P = 5 erron e i) + d f ertate A EA cos Y (A EN 
E SH ER" /E/2N 7 o JE 
O -T/2 O 
(3.68) 
2. Suare Wave 

Consider now the second special case in which v(t) 

is a periodic square wave as shown in Fig. 2.6. In Chapter II 


we showed that 
| (EN ee = dls 
2 
vl 
0 for k = even 
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so that equation (3.59) becomes 


E Ro 
; cos k21 (2352) 
a 2 





MA) = PEN. k=1 k 
g O © 
É Cam 
2=1 2° 
Zo - -a, 2 
e A + (49) 
= vV2E/N 
O ee, 
T /6 
0 2 > cs) + 6 (25%) $) (3.69) 
Therefore 
ZTA 
il | y A 1 T/2 mcos——= 
Ea SE 7 erfcx (————) T I. (mt 
| VEN, ey 
Br We ns 
erfx (—— - YZE/N (1 - 6%) + rar (3.70) 
VE/EN, T T 
Observe that equations (3.35) and (3.62) involve 
summations whose indices run from 1 to infinity. In practice, 
RES not possible to compute infinite sums. However, due to 
the fact that for all reasonable signals v(t), the magnitude 


of the Vy coefficients gets smaller as k gets bigger, this 
w anite sum can be truncated without introducing significant 
computational error. 


Fortunately for triangular and square wave signals 


considered as special cases in this chapter, the argument of 


Ee 


summation in equations (3.35) and (3.6%) turned out to ber 
such a form that it could be expressed in closed form. 

For the sine wave, summation disappeared due to the fact 

that its coefficients Vy exist only for k = +l. These results 
in the form of graphs are presented and discussed in Chapter 


IV. 
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IV. DISCUSSION OF GRAPHICAL RESULTS 


This chapter presents graphical results obtained by 
applying the analytical results of the previous chapters to 
specific examples and carrying out the required computations 
on the computer. The plots are intended to display receiver 
performance in terms of probability of error (P,) as a function 
of SNR. Some curves of Po versus SE threshold level 
have also been included in order to show the dependency of Pa 
on the threshold which in turn depends on SNR. The graphical 
results displaying performance of the optimum receivers are 
presented first, followed by graphical results displaying 


performance of the suboptimum receivers analyzed. 


Pork earHICAL RESULTS FOR OPTIMUM RECEIVERS 


l. Receivers Operating on Signals with Nonzero DC 
Component (Uniform p-d- f- on À) 


In Chapter II, using statistical communication theoretic 
principles, a simple integrator receiver (Fig. 2.2) was derived 
that could discriminate between signal plus noise and noise 
only hypotheses provided the signal had a non-zero DC component. 
Recall that the receiver derivation and the subsequent per- 
formance evaluation carried out was made possible by the 
two-term approximation on the exponential appearing in Equation 
(2.19). Equivalently, low SNR conditions were assumed so that 


the P results are only valid for low SNR's. 
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Figure 4.1 shows the performance of this receiver as 
a function of SNR. This graph was obtained through numerical 
evaluation of Equation (2.40) under the assumption that v(t) 
Has a strong DC component in comparison to its harmonics, 
Serer lor probabilities, and uniform p.d.f. on the r.v. A. Due 
to the constraint imposed by Equation (2.20), Equation (2.40) 
Memvalid up to an SNR of about -4.1 dB. Thus Fig. 4.1 is 
plotted only up to the SNR of O ae 

Observe that PL is at a very high level of about 0.3. 
However, it decreases as SNR increases, as expected. While 
it is clear that no receiver could operate with an error 
probability of 0.3, the figure does give an indication of the 
performance level that can be expected at low SNR's. For 
high SNR's, it is expected that Gë will continue to decrease 
to reasonable levels. However, for high SNR's, the receiver 


Eo. 2.2 is no longer optimum. 


2. Receiver Operating on Arbitrary Signals (Nonuniform 
Bed. ka Ole) 


The two-term approximation on the exponential appearing 
Men uation (2.19) and application of the nonuniform p.d.f. on 
A given by Equation (2.1) resulted in the receiver of Fig. 2.4 
that could operate on signals having either zero or non-zero 
DC component. Observe that as m > 0, the receiver of Fig. 
2.4 becomes equal to the receiver of Fig. 2.2. Two specific 
signals, a sine wave and square wave, both having zero DC 


component, were analyzed. 
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Figure 4.2, obtained through numerical evaluation 
of Equation (2.80), shows the performance of this receiver 
as a function of SNR when v(t) is a sine wave. Three differ- 
ent values of the independent parameter m were considered. 
Recall that m controls the level of uncertainty about the time 
delay A. As expected, as m increases, ra decreases for a 
given SNR. For example, æa SNR of -5.2 dB is required in 
order to achieve a Po OT E m= 2. With m = 90, we need 
an SNR of -9 dB to achieve the same Pa. Note that the m = 0 


0 


case has not been considered here due to the fact that m 
corresponds to a uniform p.d.f. on A The receiver for TENCIS 
case (Fig. 2.4) is unable to detect signals having zero DC 
component, resulting in Po equal to 0.5. 

Curves shown in Fig. 4.3 are the result of nunca 
evaluation of Equation (2.90) in which three different values 
of the independent parameter m havi DEC nE Omo C E cano 
m = 0 case has not been considered here either due to the 
Same reasons presented for the Sine wave case. 

Observe that Equation (2.90) involves the computation 
of infinite sums which can not be carried out in practice. 
However, due to the fact that for all finite power signals 
v(t), the magnitude of the Vi coefficients gets smaller as k 
gets bigger and thus, truncation is possible. Inclusion of 
more coefficients in Equation (2.90) corresponds to increased 
accuracy in the results. Furthermore, ae changes as k is 


incremented, however up to a point only as further increase in 
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Figure 4.3 Performance (Square Wave) 
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k beyond this limit do not result in an appreciable change 


in Po: This phenomenon is illustrated in Table I. 


TABLE I 


PROBABILITY OF ERROR VS * OF COFFICIENTS SUMMED 


( SNR = —10 DB) 


S OF COFFICIENTS i PROBABILITY OF 


SUMMED . ERROR 


0.4351 2437D+00 
O. 44433947D+00 
0.444325590+00 
0.444325430+00 
0.44432542D+00 
0.44432542D+00 
0.0 44432542D+00 
00444325420+00 
0.444 325342D+00 
0-44432542D+ 00 


OVONOUF WHR 


ec 





We see that P keeps on decreasing until the number of 
coefficients summed is 5. For k > 5, Po remains unchanged. 
Thus, for this example the computation of Equation (2.90) was 
Carried out with sums truncated to 5 terms. The performance 
curves for the square wave signal are shown in Fig. 4.3 
for low SNR conditions. With m = 2, a SNR of -4.1 dB is 
required in order to achieve a Pa Of 0.4 while for the same 


E with m = 90, the SNR required is -8.8 dB. 


ER GRAPHICAL RESULTS FOR SUBOPTIMUM RECEIVERS 
In Chapter III, two suboptimum receivers were analyzed. 
The design approach used in that chapter was different in the 


sense that receivers were derived via strictly heuristic 
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means. Therefore, the Optimum threshold setting = or wer 
receiver was not readily obtainable. However, a computational 
approach was used to solve the problem of properly setting 
the threshold prior to evaluating FE as a function of SNR: 
l. N-Correlator Receiver 

The first suboptimum receiver analyzed in Chapter III 
was the N-Correlator receiver of Fig. 3.1, and its performance 
with v(t) a periodic a pulse was evaluated. The 
graphical results corresponding to the. numerical evaluation 
of Equation (3.46) are shown in Figs. 4.4 through 4.12. As 
mentioned earlier, uSing the computer, Equation (3.46) has 
been evaluated as a function of the threshold y for fixed SNR 
values starting at -10 dB and increasing in increments of 5 
dB. For a given SNR, the threshold value for which En is 
minimum corresponds to the optimum threshold for that part cu 
lar SNR.. These minimum points corresponding to each SNR have 
been used as optimum threshold values in Equation (3.46) and 
the resulting P, has been plotted in Fig. 4.12 as Po VS SNR. 

Observe that Equation (3.46) involves the previously 
encountered parameter m and N, namely the number of correlators 
used by the receiver. Two different values of each parameter 
m = 0, 90 and N = 2, 16 have been used in the computations 
in order to evaluate their effect on receiver performance. An 
inspection of Fig. 4.4 (m = 0, N = 16) and Fig. 4.6 (m = 30, 
N = 16) reveal that for a given N, m has no significant effect 
on Ear On the other hand Fig. 4.4 (m = 0, N = 16) and Fig. 


4.8 (m = 0, N = 2) reveal a small effect of Non Po for a 
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given value of m. The actual effect of Non Po is shown in 
Fig. 4.12 where a has been plotted as a function of SNR for 
two values of N. Observe that the higher value of N vields 


better performance. To achieve a Po OE SCH 


with N = 2, the 
SNR required is ~ 13.5 dB whereas for the same Da with N= 16 
we need an SNR ~ 13.2 dB. This result is once again expected 
as using few correlators might tes in a small output on 
the correlator branches nee) received in rí(t) has a 
delay very much different than the locally generated v(t). 

While it is possible to "see" the reason for the 
independence on m in Equation (3.46), a heuristic explanation 
of this observed result follows the following argument. The 
receiver of Fig. 3.1 correlates the incoming signals with local 
replicas of v(t) delayed by multiples of T/N. Depending on 
the actual delay in v(t), upon reception of r(t), one of the 
correlator branches will produce the largest output (or there 
may be a tie between two branches). The output produced, 
however, namely which correlator branch is largest, should 
be independent of the level of uncertainty about A. Clearly, 
as m increases, it is easier to predict which correlator 
branch output will be largest. But this is inconsequential 
as the receiver is only trying to detect the presence of v(t). 

2. Estimator-Correlator Receiver 

The second suboptimum receiver analyzed in Chapter III 

was the Estimator-Correlator of Fig. 3.4. Its performance 


for v(t) either a periodic sine wave or a square wave was 
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PE VS THLD 
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Figure 4.4 Pe vs hlad (mae, N= 16), (GNR = -0 >, 
QO > ae) 
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Figure 4.5 
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considered as a function of SNR. As in the case for the N- 
correlator receiver, this receiver does not have an optimum 
threshold that is obtained from the theoretical development. 
However, a procedure similar to that used for the N-correlator 
receiver has been used to obtain the optimum threshold setting. 
The PL has been evaluated as a function of SNR except that 
here, plots showing Po as a function of y have not been included. 

Figures 4.13 through 4.16 have been obtained through 
numerical evaluation of Equations (3.67) and (3.70). Two 
curves in each figure have been plotted to show the effect 
of the parameter m on the receiver performance. Observe that 
Equations (3.67) and (3.69) involve a parameter a which can 
only take on values iTYN, where i = 1,2,...,N. Observe that 
i = N corresponds to a = T or equivalently a = 0 since v(t) 
MS periodic. 

Two values of a have been considered for each of the 
signals, namely the sine wave and square wave, in order to 
determine its effect on receiver performance. An inspection 
Bene. 4.13 and Fig. 4.14 displaying the performance of the 
receiver for v(t) a sine wave, or inspection of Fig. 4.15 and 
Fig. 4.16 displaying the performance of the receiver for v(t) 
a square wave, it can be observed that a does not have any 
significant effect on Es for either case. The logic support- 
ing this observation on receiver performance as a function of 
a follows from the fact that a is just the estimate of the 


time delay of the signal, which is obtained by choosing that 


E 


correlator branch that produces the largest Output IR 
actual value of a is irrelevant in so far as the performance 
of the receiver is concerned. The receiver simply uses the 
estimate of the time delay to produce a decision based on a 
correlation operation. However, one of the factors that might 
improve or degrade the error probability is the accuracy of 
the estimate of the time delay of the signal, which is 


itself a function of N and SNR. 
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Figure 4.14 Performance (Sine Wave), (Alpha = T/2) 


100. 


PERFORMANCE (SQUARE WAVE 
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Figure 4.15 Performance (Square Wave), (Alpha 
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VIE ONETUSTONS 


In this thesis, the analysis and design of noncoherent 
digital receivers that optimally (and in some cases subopti- 
mally) detect an arbitrarily shaped periodic waveform with 
random time delay, has been carried out. The probability 
gensity function(p.d.f.) of Ee random time delay À is assumed 
to be of some known form. A uniform p.d.f. and a (parameter 
varying) non-uniform p.d.f. have been considered and their 
effect on the receiver performance has been studied. 

The optimum (in minimum probability of error, Pa’ sense) 
receivers have been designed on the basis of likelihood ratio 
test under low SNR assumptions. The design of suboptimum 
receivers is based on heuristic approaches that intuitively 
yield structures or algorithms that yield good performance. 

The analysis carried out in this thesis, allowed for the 
possibility of unequal prior probabilities of occurrence of 
the hypotheses, however all examples and computations assumed 
equal prior probabilities for simplicity. 

The low SNR assumption allowed certain mathematical 
complexities to be overcome. Therefore graphical results 
involving the performance (Pu) of optimum receivers are valid 
and therefore presented only for SNR values less than or equal 
to -4 dB. Operating at SNR values less than or equal to 0 dB 
is the basic requirement of covert communication and may 


also be necessary in certain radar applications. 
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The optimum receiver operating on the signals that have 
non-zero DC component, detects the signal with higher level 
of error probability and it decreases with the increase of 
SNR as expected. While it is clear that no receiver could 
operate with such an undesirable level of error probability, 
the graph displaying the performance (Rag) does give an indi- 
cation of the performance level that can be expected at low 
SNR's. In the region of high SNR, it is expected that errol 
probability will continue to decrease to reasonable levels. 
However, for high SNR's this receiver is no longer optimum. 

The second optimum receiver operating on arbitrary sig- 
nals does also yield poor performance. But here again we 
must not forget that the receiver is operating in the low 
SNR region. The graphs displaying the performance of this 
receiver for the sine wave and square wave give an indication 
of the performance that can be expected at low SNR's. Also, 
they reveal an improvement in the performance with an in- | 
crease in the parameter m. It is concluded that for high 
SNR's, this receiver is also no longer optimum. 

One of the receivers designed via heuristic means was a 
(suboptimum) N-correlator receiver. Its performance was 
derived and tested for the case in which v(t) was a periodic 
triangular waveform. While this receiver is unable to detect 
a sine wave, it is capable of detecting a square wave or some 
other periodic waveform. The other receiver designed via 
heuristic means, namely the Estimator-Correlator receiver, 


was analyzed in terms of its P, as a Func ION OMS NA 
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From the graphical results of the performance of the N- 
correlator receiver, it is concluded that receiver performance 
gets better as N, the number of correlators, gets bigger. 

It is further concluded that m, the parameter that determines 
the uncertainty about the time delay à, does not have any 
effect on the receiver performance. In general, it was 
observed that except for the case of the N-correlator receiver, 
the parameter m had a EE effect on the performance 

of the receivers. 

while a number of practical problems have been addressed 
in this thesis, related problems and potential studies abound. 
For example, the receiver obtained using a three-term approxi- 
mation in the exponential of Equation (2.19) needs to be 
analyzed further in terms of performance for both uniform and 
non-uniform p.d.f. on A. Furthermore, other (suboptimum) 
receivers need to be proposed and analyzed, such as simple 


energy detection receivers. 
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APPENDIX A 


DERIVATION OF MEAN AND VARIANCE OF R.V. A 


The p.d.f. of r.v. A (time delay) given by Equation NLM 
is 


21 
m cos A 
e 


Wi, tay = 7/2], 0<m<o (A.1) 


£, (Arm) 
Substitution of Equation (2.41) into Equation (A.1) yieles 


) I. (m) cos pan 
: o ne = 
f Azm) TI, m) ASA O ES (A.2) 


so that the expected value or mean of r.v. \ can be expressed 


as 


772 
E(A) = f AL, (Am) da 
-T/2 
1 S pe 
a _ A ee cos (p21\/ Ex 
TI, (m) p==0 p -T/2 
= 0 (ATEN 
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mece. the product of À (odd function) and cos(p2mA/T) (even 
function) is an odd function and integral of an odd function 
produces zero. 


The variance of A is given by 


Do 
E{\} = f ACE, (Arm) dd 
-T/2 


122 
1 f ,2¿McOs 2774 /T 


ner dA LEON) 
EE 


since we have a zero mean random variable. 
The variance is clearly a function of m. However, the 
Mmtegral of Equation (A.3) is not readily evaluated in closed 


form. The computer has been used to evaluate 


var{\} _ la f <2 eMCOSTX gx (A.4) 


and it has been clearly demonstrated that Var{A}/T“ becomes 


smaller as m grows. 
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APPENDIX B 


INTERPRETATION OF A SMALL EXPONENTIAL 


IN EQUATION (2.22) AND EQUATION (27 6G) 


In Equation (2722) and Equation) 


that 


L 


0 A 
f Cit vice) cl 
0 


€ << IL 


a 
N 
O 


This appendix is devoted to analyzing under what 


it has been assumed 


CERA 


conditions 


this assumption valid. Observe that g is a r.v. so that 
Equation (B.1) must be analyzed in a statistical sense. Rec 
that 
E (BZ) 
since the noise has been assumed to be zero mean and 
E 
o 2 2 = 2 
EI A E (B.3) 
o 0 O 


where E is the signal energy for D Sea To: 


Thus, provided E/N, << l, e under either hypothesis on 
the average takes the value zero or a value close to it. The 
excursions of E away from the mean can be obtained by evadeg 


ing the conditional variances. Since 
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T 


2 2 S 2 
Efe |H} = Et [y= J n(t)v(t-d) dt] *} 
0 
a 
4 0 
A Ge A Ea 
N 0 
Ha 
2 2 2 
Ze v (t-))dt = — E (B.4) 
SL a No | 
and 
To "o 
EOS le. 3 - EI Li f v7 (t-d) dt + + f n(t)v(t-A) dt] 7} 
o 0 Sch 
SS 
= Ele)” + 205 fs n(t)v(t-a)dt) 
O N^ 0 
O 
To 
+ (É TI  n(t)v(t-a))*) 
o 0 
2E 2 FE 
= IN) Enz (Des) 
O O 
Thus 
2E 
Var{e|H,} = = = Varfe|H,} (B.6) 
ee IH, 
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Therefore, provided that E/N, << 1, e remains very close 
to its average value which is either zero or very close to 


zero. Note also that in a more rigorous sense, we can bound 


Die >S|H,} e ENE OS (BR 
If we set ô = l, then 
2 2E 
Ple >1|H,) Ss (B.8) 
O 
and 
Ple” >1|H,] < ey t ce (BAS) 
O O 
In both cases if 2E/N << 1, then Ple >1 or € <-1} is also 
a as = 
small. Therefore, a small E/N, guarantees that the assumption 


of a small exponential in Equations (2.22) and (2.60) is 


valid with a high degree of probability. 
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